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1. Introduction 

One of the most favoured possible extensions 
of the Standard Model (SM) is supersymmetry 
(SUSY), which renders natural the electroweak 
mass scale [1] and accommodates grand unifica- 
tion of the particle interactions [2] . If i? parity is 
conserved it also provides a promising candidate 
for astrophysical cold dark matter, which might 
be the lightest neutralino, x? [3], or the gravitino 
[3l4l5j . SUSY also predicts the appearance of a 
relatively light Higgs boson [617] . and may pro- 
vide a welcome correction to the SM prediction 
for the anomalous magnetic moment of the muon, 
(g~2)^ [SMTU]. 

However, even the minimal supcrsymmctric ex- 
tension of the SM, the MSSM [11], boasts over 
100 free parameters, mostly associated with the 
mechanism of soft SUSY-brcaking. Hence sim- 
plified scenarios with particular restrictions on 
the pattern of SUSY-breaking are often studied. 
One example is the constrained MSSM (CMSSM) 
|12|13|14|15| . in which the soft SUSY-breaking 
gaugino masses mi/2, scalar masses mo and trilin- 
ear couplings Aq are each assumed to be universal 
at the grand unification (GUT) scale, and tan /? is 
unconstrained. This leads to four effectively-free 
parameters, if the gravitino is assumed to be suf- 
ficiently heavy and/or rare that its cosmological 
decays and its mass are irrelevant. Another pos- 
sibility is to relax the universality constraint for 
common soft SUSY-breaking contributions to the 
Higgs masses, yielding the NUHMl [16] with five 
effective parameters in addition to the gravitino 
mass that is assumed to be irrelevant. 

Alternatively, additional assumptions may 
be imposed, as in minimal supcrgravity 
(mSUGRA) I17I18I19I20I2T] . in which there is a 
specific relation between the trilinear and bilinear 
soft SUSY-brcaking parameters and the universal 
scalar mass: Aq = Bq + mo, and the gravitino 
mass is set equal to the common scalar mass 
before renormalization: = 'tt-o Q- Hence 

mSUGRA has just 3 free parameters, namely 
^1/2, 'Tio and Aq, and tan/3 is now fixed by the 
radiative electroweak symmetry breaking condi- 

^See, for example, remark (b) following equation (16) 
of [19]. 



tions [22. Further, there is a restriction on mo if 
the lightest supcrsymmctric particle (LSP) is the 
lightest neutralino, Xij assume here for the 

reasons discussed below. An intermediate sce- 
nario is the very constrained MSSM (VCMSSM), 
in which again ^o = -^o + mo but the gravitino 
mass is left free, so there is no restriction on 
mo [23j . If m3/2 is sufficiently large, and/or the 
gravitino abundance is sufficiently low (as we will 
assume here), there are no related cosmological 
constraints, and the VCMSSM also has effec- 
tively 3 relevant free parameters, but they are 
less constrained than in mSUGRA. 

We have previously published frequentist anal- 
yses of the CMSSM and NUHMl parameter 
spaces |24I25I26I27I28| . implementing the exper- 
imental constraints from low-energy electroweak 
precision data and the lightest Higgs boson mass 
Mfi as well as the lower limits from the direct 
searches for SUSY particles at LEP, fitting the 
measured value of the anomalous magnetic mo- 
ment of the muon, (g — 2)^, B physics and the 
cosmological dark matter density, Q^h^, assum- 
ing that the LSP is the lightest neutralino, x^- 
In this paper we extend these analyses to include 
the VCMSSM and mSUGRA which have, as dis- 
cussed above, one or two additional constraints on 
the pattern of soft SUSY-breaking, respectively. 
An early analysis in these scenarios can be 
found in jl5j . 

In the case of the VCMSSM, in which = 
Bq + mo but the gravitino mass is free, we as- 
sume that the gravitino is sufficiently heavy or 
rare that the dark matter is composed of neu- 
tralinos and the cosmological effects of its decays 
are unimportant. Under these assumptions, as we 
show below, imposing the neutralino dark mat- 
ter constraint docs not increase substantially the 
of the global minimum, which is ^ 1.2 higher 
than in the CMSSM, but we find 68 and 95% 
confidence- level (CL) ranges of mj^/2 and mo that 
are more restrictive than those found previously 
in our analyses of the CMSSM and the NUHMl. 

However, the preferred part of the VCMSSM 
parameter space has m-^^o > mo, which within 
mSUGRA would imply that 7TI3/2 = mo < m^j^o, 
so that the lightest neutralino would be unstable, 
and the dark matter would be composed of grav- 
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itinos. In such a case, the usual calculation of the 
ncutralino dark matter density Q.^Jt? would be in- 
applicable, and one should, instead, consider the 
constraints on the decays of long-lived neutrali- 
nos or other sparticles into gravitinos that are im- 
posed by the cosmological abundances of light ele- 
ments [5|29l30j . In fact, these constraints arc suf- 
ficiently strong to exclude, within mSUGRA and 
assuming a standard cosmological evolution, the 
otherwise preferred regions of parameter space 
with Too = TO3/2 < TO^o and hence gravitino dark 
matter [29] H. Thus, the surviving region of the 
mSUGRA parameter space has TO3/2 = too > 
TO^o, corresponding to ncutralino dark matter. 

We find within the mSUGRA model two dis- 
tinct regions with local minima of the func- 
tion, each of which is significantly worse than in 
the VCMSSM, namely Ax^ - 7 or 11 with 19 dof, 
corresponding to a goodncss-of-fit of just 6.0% 
(or 2.3%) in the mSUGRA hypothesis. This may 
be contrasted with the cases of the CMSSM and 
the NUHMl studied in |24I25I26I27| . where the 
best-fit parameters are consistent with the cur- 
rent experimental constraints at the level of 32% 
(31%) fit probability, and with the VCMSSM case 
(20 dof, 31% fit probability). One of the regions 
preferred in mSUGRA has larger values of both 
TOi/2 and Too than in the VCMSSM, CMSSM and 
NUHMl, and the other has larger too but small 

TOi/2. 

2. Notations 

Before describing our analyses of the VCMSSM 
and mSUGRA in more detail, we first specify 
our notations, since different conventions for the 
MSSM superpotential couplings and the trilin- 
ear and bilinear soft SUSY-breaking terms are 
used elsewhere in the literature, including in |24j . 
Our conventions here follow those specified, e.g., 
in [32| , according to which the superpotential in- 
cludes the terms 

W 9 YeeHiLE" + YdeHiQD" (1) 
+ YueH2QU^ + HeHiH2 , 

^See 1311 for possible effects of non-standard cosmological 
histories that might invalidate these arguments against the 
gravitino LSP scenario. 



where the Yj are Yukawa couplings, e is the 
antisymmetric 2x2 tensor with £12 = -1-1, 
Hi 2, L, E'^,Q, D'^ and If^ are superfields, and /i 
is the Higgs supermultiplct mixing parameter. 
The corresponding trilinear and bilinear SUSY- 
breaking terms in the effective Lagrangian are: 

C 3 ~{AtYtth2qt^ + ...)- ^iBehih2, (2) 

where the lower-case letters denote the scalar 
components of the corresponding superfields. 
Within this convention, sin 2/3 = —2Bijl/ {m\ + 
m| -I- 2//^) at the tree level, where toi_2 is the 
soft SUSY-breaking mass of i?i,2, the left- right 
mixing term in the stop mixing matrix is to|^ = 
—mt{At + /zcot/3), and the one-loop renormal- 
ization of the trilinear coefhcient has the form 
dAt/dt 9 - f 5IM3 + QY^At + (where M3 de- 
notes the soft SUSY-breaking parameter in the 
gluino sector, and At is the trilinear Higgs-stop 
coupling.) These choices unambiguously deter- 
mine the sign conventions for A and B. 

Within this convention, ^0 = + toq before 
renormalization in mSUGRA with its minimal 
(flat) Kahler potentialH. Additionally, as already 
mentioned, the choice of a minimal Kahler po- 
tential also imposes the condition TO3/2 = TOq on 
the gravitino mass before renormalization, so that 
mSUGRA has just three independent parameters. 

3. Details of the evaluation 

Our analysis has been performed using the 
MasterCode |24l25l26l27l28| . We sample the 
VCMSSM and mSUGRA parameter spaces us- 
ing a Markov Chain Monte Carlo (MCMC) tech- 
nique similar to that used in our previous analy- 
ses of the CMSSM and NUHMl |24I25I26I27| . We 
evaluate the global likelihood using a function 
constructed by combining the likelihoods for the 
experimental constraints from electroweak preci- 
sion data, the anomalous magnetic moment of the 
muon, {g — 2)^, B physics, the astrophysical cold 

•^Note that many of the publicly available packages such 
as Sof tSusy 1331 use the opposite sign convention, as may 
bo ascertained by comparing the signs of the gauge and 
Yukawa contributions to the RGEs of the A parameters. 
In the notation of those codes, the mSUGRA boundary 
condition would bo Aq = Bq — mo . 
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dark matter density, il^h^, and searches for the 
hghtest Higgs boson and supersymmetric parti- 
cles, in exactly the same way as described pre- 
viously |24I25I26] . The most significant change 
in our numerical treatments of these observables 
since [26] is in {g — 2)^. for which we use use the 
estimate a^uSY = (28.7 ± 8.0) x 10"!° [10]. 

The numerical evaluation within the 
MasterCode combines the following theoretical 
codes. For the RGB running of the soft SUSY- 
breaking parameters, it uses SoftSUSY [33] . 
which is combined consistently with the 
codes used for the low-energy observables: 
FeynHiggs [7I34I35I35] is used for the evalua- 
tion of the Higgs masses and a^^^^ (see also 
[37138139140] ). For flavour- related observables we 
use SuFla |41|42| and Superlso i43j44y, and for 
the electroweak precision data we have included 
a code based on [45l46j . Finally, for dark-matter- 
related observables, we use MicrOMEGAs |47i48i4g] 
and DarkSUSY |50I51| . We make extensive use 
of the SUSY Les Houches Accord |52l53j in the 
combination of the various codes within the 
MasterCode. 

Our MCMC samphng of the VCMSSM param- 
eter space comprises some 30,000,000 points. The 
neutralino CDM constraint on fi^/i^ [54] and the 
Higgs mass constraint |55|56j were applied after 
the sampling, allowing the effects of these two 
constraints to be studied separately. In the case of 
mSUGRA, about 17,000,000 of the MCMC points 
from the VCMSSM sample survive the mSUGRA 
constraint mo = 7713/2 > "Z-j^o, and we again ap- 
plied the neutralino CDM constraint on il^h^ and 
the Higgs mass constraint a posteriori. 

4. Analysis of parameter planes 

We start our analysis with the results in the 
(mo,mi/2) planes. Fig. [1] displays the global like- 
lihood functions in the VCMSSM (left panels) 
and mSUGRA (right panels). In each case, the 
upper panel shows results before the fty^h'^ con- 
straint is applied, and the lower panel displays 
the effects of imposing the fl^h^ constraint |j. In 

*We recall that the constraints due to the late decays of 
massive metastable particles 1291 are taken into account 
implicitly as described above, i.e., in the VCMSSM by 



all panels, we display the points with the minimal 
values of (green stars) as well as the 68 and 
95% CL contours (red and blue), corresponding 
to Ax^ = 2.28 and 5.99. Other contours of Ax^ 
are indicated in shades of grey. 

In the upper left panel for the VCMSSM be- 
fore applying the fl^^h'^ constraint, the triangu- 
lar region at small mo and large mi/2 is ex- 
cluded because there the LSP would be charged, 
and a band extending to large mo at low mj^/2 
is excluded by the LEP Higgs constraint. The 
best-fit point is at (mo,mi/2) = (30,310) GeV, 
the 68 and 95% CL contours enclose regions of 
the (mo,m]^/2) planes that are similar to those 
favoured in the CMSSM and NUHMl [HHH], and 
there are no preferred points in the focus-point 
region at large mo. Specifically, we find that 
the focus-point region at large mo ^ 2000 GeV 
with mi/2 ~ 200 GeV is subject to a penalty of 
Ax^ - 10. 

In the upper right panel for mSUGRA with- 
out the fl^h^ constraint, we see a similar pat- 
tern, but with a larger excluded triangular region 
at large mi/2 and small mo, as the allowed part of 
the {mo, nil ^2) plane is now restricted to the re- 
gion where m3/2 = > ^n^j-o ~ 0.4mi/2- Since 
the best-fit VCMSSM point seen in the upper left 
panel of Fig. [T] lies in the region that is disallowed 
in mSUGRA, there is a new mSUGRA best-fit 
point on the boundary of the allowed region, 
with (mo,mi/2) = (110, 280) GeV. The mini- 
mum value of x^ is higher than in the VCMSSM 
model by ^ 1.5, and the Ax^ values in other re- 
gions of the (mo,mi/2) plane are correspondingly 
reduced, leading to the emergence of 'archipelago' 
of points at (m,o,mi/2) - (700,1000) GeV that 
are now allowed at the 95% CL. The fact that 
the mSUGRA best-fit point lies on the bound- 
ary of the allowed parameter space indicates that 
the restrictions in this model are disfavoured by 
current experimental data. 

We recall that it was shown in [5S] that in 
the CMSSM the fl^h'^ constraint [54] has a rel- 
atively modest impact on the preferred ranges 

assuming that the gravitino mass is high and/or its pri- 
mordial density is low, and by accepting that the de- 
cay constraints are so severe in mSUGRA as to forbid 
™3/2 = ™0 < ^j^o ™ that case. 
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Figure 1. The {mo,mi/2) planes in the VCMSSM (left panels) and mSUGRA (right panels), without 
(upper) and with (lower panels) the fl^h^ constraint \54^ , showing in each case the best-fit points (green 
stars) and the 68 and 95% CL contours (red and black, respectively). The open green star in the lower 
right panel denotes the secondary minimum discussed in the text. 



in the {mo,mi/2) plane. However, imposing the 
constraint has a dramatic effect on the 
VCMSSM fit, as we see in the lower left panel 
of Fig. □ The region allowed at the 95% CL 
is reduced to a narrow 'WMAP strip' terminat- 
ing at (mo,mi/2) - (250,700) GcV. We recall 
that similar WMAP coannihilation strips appear 
in the CMSSM for fixed values of tan/3 and Aq, 
but move across the (mo, m.1/2) plane as tan (3 and 
Aq are varied, which was why the WMAP strip 
structure was invisible in the global likelihood fit 
to the CMSSM [25]. On the other hand, we recall 
that, in the VCMSSM, tan/3 is fixed as a function 
of 77io,mi/2 and Aq, and the dependence on Aq 
is not very strong. As a result of the loss of the 



freedom to vary tan /3 independently, the WMAP 
strip structure is resurrected in the VCMSSM. 

The best-fit point in the VCMSSM with 
the fl^h^ constraint has (mo,mi/2,^o) ~ 
(60,300,30) GeV and tan/3 ^ 9. However, as in 
the previous CMSSM and NUHMl cases, rather 
larger values of mi/2, and hence m^o,mg and 
other sparticle masses, are allowed at the 95% CL. 
This VCMSSM fit has a very good value of xVdof 
= 22.5/20 (31% probability), similar to the be- 
fore applying the fl^h^ constraint, demonstrating 
that there is no significant tension between this 
and other constraints. The increase in due 
to inclusion of the LEP Mh constraint [55156] is 
^ 1.1, demonstrating that there is also no sig- 
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Figure 2. The (tan/3, planes in the VCMSSM (left) and mSUGRA (right), without (upper) and 

with (lower panels) the i^x^^ constraint f5^| /, showing in each case the best-fit points and the 68 and 95% 
CL contours. 



nificant tension between the Mh and other con- 
straints. Moreover, the fact that the fit probabil- 
ity in the VCMSSM is about the same as the value 
of 32% (x^ = 21.3/19 dof) found in the CMSSM 
indicates that applying the extra VCMSSM con- 
straint = -Bo + nii'O is certainly not a source of 
significant tension in the fit. 

Wc note that there is no focus-point region vis- 
ible in the VCMSSM when the U.-^h'^ constraint is 
applied, since it is not compatible with the tan/3 
constraint imposed by the initial conditions. On 
the other hand, there is a very narrow strip at 
mi/2 ~ 130 GeV extending to large mo, where 
the relic density is brought into the WMAP range 
by rapid annihilation through the direct-channel 
light Higgs pole. However, this strip has Ax^ > 9. 



Turning finally to the lower right panel of Fig.[T] 
for mSUGRA with the constraint applied, 

we see an evolution of the picture. Much of 
the VCMSSM 'WMAP strip' has disappeared, as 
only a vestige of it has m-^o < mo = rn^i2- Since 
the minimum value of in the VCMSSM was 
located in the forbidden region with m-^o > mo = 
m3/2, the minimum value of x^ in the mSUGRA 
region is significantly higher, specifically x^ ~ 29, 
a price Ax^ ^ 7. The best mSUGRA fit is again 
on the boundary of the allowed region. Moreover, 
as indicated by the open green star in the lower 
right panel of Fig. [TJ a second local minimum with 
X^ ^ 33 (which is therefore allowed at the 95% 
CL) can be found along the light Higgs rapid- 
annihilation strip with mi/2 ^ 130 GeV and 
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900 GeV < mo < 2500 GeV. Along this strip, 
the function is relatively insensitive to mo, 
thanks to approximate compensation between the 
contributions from BR(6 — )■ 57), {g — 2)^ and the 
forward-backward asymmetry of b quarks mea- 
sured at LEP, decreasing slightly as mo increases 
up to mo < 2100 GeV. 

Fig. [2] displays the (tan/3, mj^/2) planes in the 
VCMSSM (left) and mSUGRA (right, where the 
cut m-^o < mo was applied), both without (up- 
per) and with (lower) the Q,^h^ constraint. In the 
upper left panel, we see that a range of tan (3 ^ Q 
to 12 — similar to that favoured in the CMSSM 
and NUHMl — is preferred in the VCMSSM at 
the 68% CL, but with a best-fit value tan /3 ~ 9 
that is rather smaller. A range of slightly larger 
tan /3 - 6 to 13 is allowed in mSUGRA at the 68% 
CL, with a best-fit value of tan/? ^ 8 before im- 
posing the Vl^h? constraint. The 95% CL ranges 
of tan/3 extend to - 30, 35 in the VCMSSM and 
mSUGRA, respectively. When the ^^h^ con- 
straint is imposed on the VCMSSM (lower left), 
the ranges of tan /3 favoured at the 68 and 95% CL 
are little changed, with tan/3 ^ 9 at the best-fit 
point. We again see at higher Ax^ the rapid- 
annihilation Higgs funnel at mi/2 ^ 130 GeV, 
separated from the favoured coannihilation re- 
gion. However, when the fl^h^ constraint is im- 
posed on mSUGRA (lower right), the coannihi- 
lation region shrinks to a vestigial region with 
tan/3 30 and relatively high x^, as previously 
remarked, and the other minimum along the light 
Higgs funnel has Ax^ ~ 4 and tan/3 ~ 5 to 14, 
and also contains an area allowed at the 95% CL. 

Fig. [3] displays the corresponding (mo, tan /3) 
planes in the VCMSSM (left) and mSUGRA 
(right). In the absence of the fl^h^ constraint 
(upper panels), we see regions allowed at the 
95% CL that are restricted to mo ^ 300 GeV 
in the VCMSSM case (left) and - 500 GeV 
(~' 1000 GeV including the 'archipelago') in the 
mSUGRA model (right), the larger range being 
expected from the restriction mo > m^o. Impos- 
ing the ri^/i^ constraint on the VCMSSM (lower 
left), we find a band at low mo that can be iden- 
tified with the WMAP-compatible coannihilation 
strip in the corresponding panel of Fig. [TJ Turn- 



ing to the case of mSUGRA with the fi^/i^ con- 
straint (lower right), only a vestige of the low-mo 
band remains at relatively large tan/3, and the 
X^ in this region is '--^ 4 smaller than in the rapid- 
annihilation region at larger mo and smaller tan /3 
values H. 

Finally, Fig. [4] displays the (^o/mo, tan/3) 
planes for the VCMSSM (left) and mSUGRA 
(right). In the VCMSSM, a wide range of val- 
ues extending from Aq/itio ^ +2 down to ~ —3 
(corresponding to the case of very small mo) lies 
within the 95% CL contour, whereas in mSUGRA 
only positive values are favoured at the 95% 
CL before the fl^h^ constraint is applied. On 
the other hand, when it is applied, large neg- 
ative values of Ao/mo are disfavoured also in 
the VCMSSM, and values of Ao/mo ~ 2 with 
tan/3 ~ 30 are allowed at the 95% CL. Two dis- 
tinct populations of preferred points arc appar- 
ent in mSUGRA when the fl^h^ constraint is ap- 
plied. There is a vestigial coannihilation region 
with Ao/mo ^ 2 and tan/3 ~ 28, and the rapid- 
annihilation funnel region has Ao/mo ^ 0.4 and 
tan/3 7. 

We recall that the simplest Polonyi model of 
SUSY breaking |T7] predicts that |^o/mo| = 
3 — -\/3, a possibility that is quite consistent with 
the VCMSSM both before and after applying the 
fl^h^ constraint (left panels), but only marginally 
consistent with the 95% CL region of mSUGRA. 

5. Comparison of Likelihood Analyses 

We gather in Table [1] some important aspects 
of the best-fit points in mSUGRA, the VCMSSM, 
the CMSSM and the NUHMl (the latter being 
adapted from |25I26| . with the inclusion of the 
updated values of {g — 2)^ and m^: we find that 
the values of the different parameters at the best 
fit points in the CMSSM and the NUHMl, as 
well as the fit probabilities, remain essentially un- 
changed with respect to the analysis in [25126] ). 

^Because the rapid-annihilation region is very narrow in 
mx/2i S'Dd quite broad in mo and Aq, it is difficult to 
sample fully using the MCMC approach, even with the 
high statistics of our full sample. Moreover, the MCMC 
approach samples the input variables mi/2)"^0 and Ag, 
contributing also to the uneven sampling in the derived 
quantity tan /3 seen in the lower panels of Fig. |3] 
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Figure 3. The (mp, tan/3) planes in the VCMSSM (left) and mSUGRA (right), without (upper) and with 
(lower panels) the i^x^^ constraint f5^| /, showing in each case the best-fit points and the 68 and 95% CL 
contours. 



In the case of mSUGRA, we list the properties 
of two local minima of the function: one is 
in the coannihilation region and one in the light 
Higgs rapid-annihilation funnel with Ax^ ~ 4, as 
discussed earlier. We see in the first column that 
the mSUGRA fits have substantially higher 
than the other models, which is reflected in the 
second column by a significantly lower probabil- 
ity. The CMSSM and NUHMI gave comparable 
fit probabilities of 32 (3I)% and the VCMSSM fit 
has a probability of 31%, while mSUGRA, with 
a probability of 6.0% in the coannihilation region 
and 2.3% in the light Higgs funnel region, pro- 
vides a worse description of the data considered 
in this analysis. As a result, the mSUGRA sce- 
nario is somewhat disfavoured compared to the 



other SUSY scenarios we consider. The source 
of tension in mSUGRA is seen by comparing the 
third and fourth columns. The best VCMSSM fit 
has very similar values of mi/2 and toq to those 
in the CMSSM, and in both cases toq ^ TOi/2- 
The conflict between this preference and the cos- 
mological requirement on mSUGRA that = 
mo > m^o leads to a best fit in the coannhila- 
tion region with larger values of toi/2 and tan/3 
as well as a larger value of mg, and the other lo- 
cal minimum with small mi/2 much larger 
mo H. As seen in the sixth column, the values 



®On the other hand, the best NUHMI fit has a value of 
mi/2 similar to those in the VCMSSM and CMSSM but 
a somewhat larger value of mo , which is possible because 
an acceptable value of Q,^h^ may be found along a funnel 
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Figure 4. The (Ao/mg, tan/3) planes in the VCMSSM (left) and mSUGRA (right), without (upper) and 
with (lower panels) the i^x^^ constraint \54^ , showing the best-fit points and the 68 and 95% CL contours 
in each case. 



of tan/3 favoured in the VCMSSM, CMSSM and 
NUHMl are much smaher than that favoured at 
the best-fit point of mSUGRA in the coannihila- 
tion region. The last column shows the best-fit 
values of M?i, not taking into account the LEP 
and Tevatron limits. It can be seen that the best- 
fit values of Mh in mSUGRA in the coannihilation 
region, the VCMSSM and the CMSSM are sim- 
ilar, and somewhat below the LEP lower limit, 
leading to increases in of 3.9 (1.1) (1-4) when 
the LEP constraint is applied. On the other hand, 
we find rather higher best-fit values of Mh in the 



in the (mo,miy2) plane at lower values of tan/3, due to 
XiX? annihilation througli direct-channel H, A poles. This 
option is not available in the CMSSM except at large tan j3 
and and it is absent completely in mSUGRA. 



NUHMl and mSUGRA in the funnel region. Fi- 
nally, we note that the largest variations between 
the different models occur for Aq (fifth column), 
reflecting the relative insensitivity of our fits to 
this parameter. 

In Fig. [S] we display the spectra in the 
VCMSSM (top) and mSUGRA (middle and bot- 
tom) with the fi^/i^ constraint applied, comple- 
menting the CMSSM and NUHMl spectra shown 
in Fig. 3 of [26]. In the case of the VCMSSM, the 
spectrum is qualitatively similar to those in the 
CMSSM and NUHMl [2§l26l . The two mSUGRA 
spectra are in the coannihilation (middle) and 
funnel region (bottom) , reflecting the coexistence 
of two qualitatively different (near-)best-fit points 
with relatively similar x^. The spectra in these 
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Model 


Minimum 


Probability 


"M/2 


mo 


^0 


tan /3 


Mh (no LEP) 


mSUGRA 


29.4 


6.0% 


550 


230 


430 


28 


107.7 




33.2 


2.3% 


130 


2110 


980 


7 


116.9 


VCMSSM 


22.5 


31% 


300 


60 


30 


9 


109.3 


CMSSM 


21.3 


32% 


320 


60 


-160 


11 


107.9 


NUHMl 


19.3 


31% 


260 


100 


1010 


8 


119.5 



Table 1 

Comparison of the best fits found in this paper within the mSUGRA and VCMSSM frameworks with 
previous results \ 251 2611 *^ CMSSM and NUHMl frameworks. In addition to the minimum value of 
in each scenario, we include the values o/ mi/2, wq, and tan/3 at all the best-fit points, as well 
as Mh (for the latter the direct bounds from LEP and the Tevatron are not included). All masses are 
in GeV units. We list two very different mSUCRA fit results with similar values: the first is in 
the coannihilation region, and the second is in the light Higgs funnel region. Note that we use here the 
convention for Aq described in the text, which differs from that in \24^ . 



regions are significantly different from each other 
and from the VCMSSM, CMSSM and NUHMl. 
This is because the coannihilation region has 
mi/2 significantly larger than in the other mod- 
els, whereas the funnel region has a significantly 
smaller and very well-defined value of mi/2 and 
relatively large values of mo. This bimodality af- 
fects directly the preferred values of mjj.o and rUg, 
and affects the other sparticle masses via renor- 
malization effects. These spectra show that the 
coloured particles are within the reach of the LHC 
for the VCMSSM and mSUGRA in the coannihil- 
iation region, whereas more integrated luminosity 
would be necessary for mSUGRA in the funnel re- 
gion (except for gluino production). In each sce- 
nario some SUSY particles should be accessible 
at an e^e~ collider, even with a center-of-mass 
energy as low as 500 GeV. 

We display in Fig. |6] the one-dimensional 
likelihood functions for (top panels) rUg and w^o, 
(middle panels) the mass differences m^^, — mg 
and mf J — m^o , and (bottom panels) BR(i3s ^ 
/i+jU^) and (Tp^, the spin-independent neutralino 
scattering cross section. The NUHMl curves are 
shown as purple dotted lines, the CMSSM curves 
as green dash-dotted hues, the VCMSSM curves 
as red dashed lines, and the mSUGRA curves as 
blue solid lines. In each plot, we display the Ax^ 
contribution of each model relative to the best-fit 
point in that model. Thus the secondary min- 



imum of the function for mSUGRA has the 
appropriate Ax^ relative to the absolute mini- 
mum in that model. 

We see that the likelihood functions for rUg 
and m^o in the VCMSSM are similar to those 
in the CMSSM and NUHMl [25l26] . with the 
most likely range of rUg ~ 700 to 800 GeV and 
m^a ^ 120 GeV. The corresponding likelihood 
functions in mSUGRA, on the other hand, are 
very different, reflecting once more the bimodal- 
ity in mi/2 and mp. In the VCMSSM, as in 
the CMSSM and NUHMl, the preferred range 
of rUg (top left panel) suggests that there may 
be good prospects for observing first hints of 
SUSY at the LHC in 2011/12, although values 
of rUg > 2000 GeV are permitted in the NUHMl 
with Ax^ ~ 5. In the rapid-annihilation strip of 
mSUGRA, rUg - 400 GeV and m^o - 55 GeV, 
putting the discovery of the gluino at the LHC 
within the reach of the LHC in 2011/12, and 
making nig potentially a powerful diagnostic tool 
for mSUGRA. However, discovery of the gluino 
would be delayed in the (more probable, but 
still disfavoured scenario compared to VCMSSM, 
CMSSM and NUHMl) coannihilation region of 
mSUGRA, where the best fit has a gluino mass 
around 1200 to 1400 GeV. 

The bimodality of the mSUGRA likelihood 
function is reflected again in the predictions for 
m^^ — mg (middle left) and nif^ — m^o (middle 
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Figure 5. Spectra in the VCMSSM (top), and mSUGRA in the coannihilation region (middle) and the 
funnel region (bottom), implementing all the constraints including that on ft^h^ . The horizontal solid 
lines indicate the best-fit values, the vertical solid lines are the 68% CL ranges, and the vertical dashed 
lines are the 95% CL ranges for the indicated mass parameters. 
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Figure 6. The likelihood functions for (top) rrig andm^o, (middle) mqj^—rrig and ruf-^ ~ m^o , and (bottom) 
BR(Bs fj,~^fj,~) (normalized to the SM estimate) and CTp^. In each case, we compare the predictions in 
mSUGRA (blue solid line), the VCMSSM (red dashed line), the CMSSM (green dash-dotted line) and the 
NUHMl (purple dotted line). 
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right) in Fig. [51 The two different mSUGRA min- 
ima yield different signs for m^^ — m^, whereas 
rrig > rriqj^ is favoured in the VCMSSM, as 
was shown previously to be favoured in the 
CMSSM |25I26| . and to a rather lesser extent in 
the NUHMl j25l26j . Thus, the sign and magni- 
tude of — nig are potential diagnostic tools 
for discriminating between different models^. We 
note also the different predictions for ^ — m^o : 
the VCMSSM predicts a very small mass differ- 
ence, as was shown previously in the CMSSM, 
and so docs mSUGRA in the coannihilation re- 
gion. However, in the rapid-annihilation region 
mSUGRA predicts mass differences that may be 
large, as was previously shown to be possible (to a 
lesser extent) in the NUHMl [26]. Thus discovery 
of a light gluino at the LHC and /or a light neu- 
tralino LSP would not necessarily imply that the 
lightest slepton would have a mass close to that of 
the LSP, and the lighter stau could be too heavy 
for an e+e~ collider with ^/s = 1 TeV. However, 
in this scenario the whole chargino and neutralino 
spectrum would be accessible at an e~^e~ collider 
with = 1 TeV, see Fig. M 

We see in the bottom left panel of Fig. [H] that 
in the VCMSSM a value of BR(B, /i+A*") 
slightly lower than the SM value is favoured, al- 
though larger values cannot be excluded. On the 
other hand; in mSUGRA in the coannihilation 
region, with its relatively large values of tan/3, 
we find a preferred value of BTi{Bs — ?> jj.^^^) 
that is slightly larger than in the SM, whereas 
a range around the SM level is favoured in the 
light Higgs funnel region. As already seen in [26] , 
the minimum of in the NUHMl is at a some- 
what higher value of BR(i?s fi^fi') than in 
the SM, whereas the best-fit value in the CMSSM 
is again slightly smaller, albeit with considerable 
uncertainty. The values below the SM predic- 
tion arise from chargino-induced Zbs penguin di- 
agrams, that yield analogous but smaller reduc- 
tions in BK{K TTvv). However, the SM pre- 



'^We recall that dijet + missing energy events due to the 
decay chain g — > q[i+q, qn — >■ q+X? a^r^ in general favoured 
when rriqj^ < mj, whereas there are expected to be a 
larger fraction of four-jet + missing energy events when 
> rUg, and different decay chains occur via q^, ii^2 
and 6i_2- 



diction has an uncertainty of about 10%, and 
0(100) Bg decays would be needed to 

match this error, so the differences between the 
mSUGRA, VCMSSM and best-fit CMSSM pre- 
dictions are probably unobservable. 

Turning to the likelihood functions for (7^^ in 
the bottom right panel of Fig. [6] (calculated as- 
suming a TT-N scattering a term Eat = 64 MeV: 
see, e.g., [57] for a discussion of the implications 
of modifying this assumption) we see that a range 
between lO^"'^ and lO^"'^ cm^ is favoured in the 
VCMSSM, whereas in mSUGRA a range between 
2 X 10~*^ cm^ and 5 x 10~^^ cm^ is favoured. 

Finally, in Fig. [7] we display the one- 
dimensional likelihood functions for in the 
VCMSSM (left) and mSUGRA (right), not in- 
cluding the direct limits from LEP and the Teva- 
trond. For each model we display the likelihood 
functions, including the theoretical uncertainties 
(red bands), which we take to be 1.5 GeV in all 
models. For comparison, we also show the mass 
range excluded for a SM-like Higgs boson (yel- 
low shading) obtained at LEP |55|56| . This limit 
is valid since the VCMSSM and mSUGRA are 
sub-spaces of the more general CMSSM parame- 
ter space, where the LEP limits have been shown 
to be valid [58159] . Values somewhat below the 
LEP exclusion are favoured in the coannihilation 
region of mSUGRA, in the VCMSSM and the 
CMSSM, entailing prices of 3.9 (1.1) (1.4), as 
discussed previously. In the mSUGRA coanni- 
hilation case, the global minimum of found 
when the LEP constraint is disregarded is in an 
isolated region at low (mg, mi/2)- When the LEP 
Higgs constraint is applied, this region is strongly 
disfavoured, and the global minimum moves to 
the green star shown in Fig. [1] located at much 
larger (mo,TOi/2) and with Ax^ — 3.9. This 
other minimum is reflected in the blue point at 
Mh = 121.1 GeV with Ax^ = 3.9 and the cor- 
responding horizontal red line visible at the top 
of the right panel in Fig. [7] Comparably large 
values of Mh ~ 117, 120 GeV are favoured in the 
funnel region of mSUGRA and in the NUHMl, 
so these models are naturally consistent with the 



*See I24I25I26I for the corresponding plots in the CMSSM 
and NUHMl. 
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Figure 7. The functions for in the VCMSSM (left) and mSUGRA (right), including the theoretical 
uncertainties (red bands), as well as the SM prediction for Mh based on a precision electroweak fit (blue 
band). Also shown is the mass range for a SM-like Higgs boson excluded by LEP (yellow shading), and 
the mass range that is theoretically inaccessible in TeV-scale SUSY (beige shading). Note in right panel 
the secondary minimum at Mh = 121.1 GeV, A^^ ~ 3.9 and the accompanying narrow horizontal red 
band. 




LEP bound on a SM-like Higgs boson. 

6. Conclusions 

We have completed in this paper the frequen- 
tist analysis of a nested sequence of variants of 
the MSSM: NUHMl 9 CMSSM 3 VCMSSM 3 
mSUGRA, discussing in each case the best-fit 
point, the minimum x^/dof, the 68 and 95% CL 
regions and aspects of the favoured ranges of par- 
ticle masses and other observables. We found pre- 
viously that the restriction from the NUHMl to 
the CMSSM does not change drastically the po- 
sition of the best-fit point or the favoured ranges 
of parameters and observables. This reflects the 
fact that the present data do not constrain signifi- 
cantly the heavier Higgs bosons of the MSSM, and 
so there is no significant tension in the CMSSM 
fit arising from that sector. Likewise, we found 
in the present paper that the restriction from the 
CMSSM to the VCMSSM does not have a large 
impact on the position of the best-fit point or on 
the favoured ranges. This reflects the fact that 



the present data also do not constrain signifi- 
cantly the trilinear soft SUSY-breaking param- 
eter Aq, and so there is no significant price 
to be paid when setting Aq = Bq + mo as in 
the VCMSSM. Fits in all three of the NUHMl, 
CMSSM and VCMSSM frameworks have good 
absolute probabilities in our frcquentist analyses. 

On the other hand, in mSUGRA two almost 
equally good best-fit points coexist, with rather 
different MSSM parameter values, significantly 
higher values of and lower absolute probabil- 
ity, somewhat disfavouring this scenario. This re- 
fiects the fact that the neutralino LSP constraint 
ruj^o < 7713/2 excludes the best fit found in the 
VCMSSM and forces instead either rather larger 
values of mQ,mi/2 and tan/3 or points along the 
narrow light Higgs rapid-annihilation funnel with 
small mi/2 but large tuq. 

In each of the NUHMl, CMSSM and VCMSSM 
there is a significant chance of observing first 
hints of SUSY in the 2011/12 run of the 
LHC, whereas this may be more problematic in 
mSUGRA, depending on whether its parame- 
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tcrs lie in the rapid-annihilation funnel or in the 
higher-mass coannihilation region. Correspond- 
ingly, early hints of SUSY at the LHC in this first 
run might further favour the NUHMl, CMSSM 
and VCMSSM frameworks over the higher-mass 
mSUGRA option in the coannihilation region. In 
the rapid-annihilation funnel region of mSUGRA 
one would expect abundant production of gluinos, 
whereas the first appearances of squarks, having 
masses in excess of 1 TeV, would only happen at 
a somewhat later stage. If no sign of SUSY parti- 
cles shows up in the early LHC run, in particular 
if the integrated luminosity (and energy) goes sig- 
nificantly beyond the currently foreseen 1/fb at 
7 TeV; the scenarios predicting a relatively low 
SUSY scale could soon come under pressure. At 
the time of writing, the ATLAS and CMS col- 
laborations are each examining some 35/pb of 
analysable data. 

Note added in proof 

Since the submission of this paper, the first re- 
sults of searches for supersymmetry have been 
published by the CMS and ATLAS Collabora- 
tions |60I61| . These change somewhat the param- 
eters of the best-fit points and the 68% and 95% 
CL regions in the models studied here, but do not 
disfavour any of them strongly j62j . 
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